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The three-dimensional XY model is studied at finite chemical potential using complex Langevin 
dynamics. An adaptive stepsize algorithm is implemented to cure the problem of runaway solu- 
tions that appears when using a constant stepsize. The validity of complex Langevin dynamics 
is tested against calculations using imaginary chemical potential and the world line formalism. 
While complex Langevin dynamics is found to work correctly at larger j3, it fails for smaller j3, 
in the region of the phase diagram corresponding to the disordered phase. Diagnostic tests are 
developed to identify behaviour symptomatic of incorrect convergence. These indicate that the 
failure does not depend on the severeness of the sign problem, but has a different origin. 
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1. Introduction 

Field theories with a complex action are difficult to treat nonperturbatively because the weight 
e~ s = \e~ s \e ,(p in the partition function is not real. Standard numerical approaches based on a 
probability interpretation and importance sampling will then typically break down, which is com- 
monly referred to as the sign problem. This is a particularly pressing issue with regards to the 
determination of the phase diagram of QCD in the plane of temperature and chemical potential. A 
comprehensive review can be found in Ref. [jl|] 

Complex Langevin dynamics offers the possibility of a general solution to this problem [§, ||, 
In this formulation, fields (j) are supplemented with a fictional time-like dimension, #, and the 
system evolves according to the stochastic equation 
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In the case of a complex action, the fields are complexified as — > R + i<p l , and the Langevin 
equations read (using real noise) 
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where the noise is Gaussian, (Tfc(#)) = 0, = 2<5xy<5(# - ■&'). In the limit that 

# oo noise averages should become equal to quantum expectation values. Since the action is 
used to compute the drift terms, but not for any importance sampling, complex Langevin dynamics 
can potentially avoid the sign problem. 



2. XY model 

Motivated by previous studies of complex Langevin dynamics for QCD in the heavy dense 
limit [||], we consider the XY model at finite chemical potential [ph. This theory is closely related 
to the Bose gas, for which complex Langevin dynamics is known to work well [f7], [Sp _ The XY 
model at finite chemical potential has the action 

2 

S = -J3£ £ 008(0, - 0. r+ v - ifiSyjo), (2.1) 

x v=0 

where < ty x < 2n. The theory is defined on a lattice of volume CI = N t Nj, with periodic boundary 
conditions. The chemical potential is coupled to the Noether charge associated with the global 
symmetry, <p x — > (j> x + a, and is introduced in the standard way [S]. The action satisfies S*(pl) = 
S(—jJ.*) and at vanishing chemical potential the theory is known to undergo a phase transition at 



p c = 0.45421 [10] between a disordered phase when j8 < j8 c and an ordered phase when j3 > j8 c . 
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The drift terms in the complex Langevin equations are given by 

K x = "J3 £ [sin(0 x R - cosh(0i - - jU5 V) o) (2.2a) 

V 

+ sin(0 v R - cosh($ - A ! _v + MM] . 
^ = -i 3 1 t cos (^ " <k%) sinh (^ " - MM (2.2b) 

V 

+ cos(0* - sinh(^ - ^ + At5 v ,o)] • 

By choosing an imaginary chemical potential /J. = ijij, the action becomes purely real and standard 
algorithms can be applied (here real Langevin dynamics is used). The behaviour at /I 2 > can then 
be assessed by continuation of the behaviour found at \i 2 < 0. 

3. Adaptive stepsize 

In order to integrate the equations, the Langevin time needs to be discretized as # = en. When 
real Langevin dynamics can be employed (e.g. when jx 2 < 0), the drift terms are bounded and a 
fixed stepsize is sufficient. However, in the case of complex Langevin dynamics, the drift terms are 



unbounded (see Eq. |2.2| ) and numerical instabilities are encountered when the forces become large. 
These cause the system to diverge along "runaway" solutions. 

This can be cured by adjusting the stepsize when the configuration approaches a divergent 



trajectory [Hip. The adaptive stepsize is implemented by monitoring the maximal force term, 

K™ x = max \K?(n) + iKl(n)\ , 

X 

and at each update defining the stepsize to be 

e n = mini e,e^— — ->. (3.1) 

Here e is the desired target stepsize and (K max ) is either precomputed or computed during the 
thermalisation period. 

All observables are analyzed over equal periods of Langevin time and weighted with the step- 
size to ensure correct statistical significance, 

{0) = l^n t (32) 

with the total number of updates such that e n is constant. 

While with a fixed stepsize it is practically impossible to generate a thermalised configuration, 
we found that with an adaptive stepsize instabilities are completely eliminated. The same result 



applies to heavy dense QCD [|11|]. 

4. World line formulation 

A useful feature of the XY model is that it can be rewritten exactly without a sign problem in 
terms of world lines. This dual formulation can be efficiently simulated using a worm algorithm 
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[10]. The partition function can be expanded into a sum over Bessel functions using the identity 

CO 

ft cos 4 



M 



k=-o 



which allows the partition function to be written as 

[k] x,V V v / 



(4.1) 



(4.2) 



The sum over [k] indicates a sum over all possible world line configurations. The action, (S) = 
— j6 ^jj^ , can be computed from 

<*>^[w-t]) w ,' 

where the average is taken over world line configurations. 
5. Comparison 

To assess the validity of the results from complex Langevin simulations, a good starting point 
is to compare with results from imaginary chemical potential simulations at \i 2 < 0. In Fig. |l] 
it can be seen that at large coupling j8 = 0.7 the action density is continuous over the boundary 
at n 2 = between real and complex Langevin dynamics. The results at j8 = 0.3 are in contrast 
to this. Here, the action density between the two regions are not in agreement, exemplified by 
the fact that at jx 2 = the two simulations give different results for (S) /&. These observations are 
corroborated by the results from the world line formalism, which are consistent with the predictions 
from imaginary chemical potential calculations. 
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(a) j8 = 0.3 



(b) = 0.7 



Figure 1: Real part of the action density around /I 2 ~ at low and high j3, with the world line result for 
comparison. The lattice volume is 8 3 . Note that with an imaginary chemical potential the theory is periodic 
under (j) — > (j> +2n/N T which yields a Roberge-Weiss transition at ji[ = n/N T , similar to what is found in 
QCD. 

In order to investigate this further, we have studied a large number of parameter values in the 
/3-jU plane. The disagreement between complex Langevin dynamics (cl) and world line formal- 
ism (wl) can be quantified by the relative difference between the expectation values of the action, 
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according to 



AS: 



(<S)wi — (5)ci 



(S) 



(5.1) 



wl 



This is plotted in Fig. || in the /3-/I plane, with the phase boundary taken from Ref. [|I0|]. It can 
be seen that the breakdown of complex Langevin dynamics is highly correlated with the phase 
boundary of the theory: complex Langevin dynamics is in agreement in the ordered phase but 
begins to break down at the boundary and is in complete disagreement in the disordered phase. 




0.5 1 



Figure 2: Colour plot of relative difference AS between action density as computed using complex Langevin 
dynamics and world line formulation. Also shown is the phase boundary /3 £ .(jli) between ordered phase (large 
j3) and disordered phase (small j3). 



6. Diagnostics 

To identify the manner in which complex Langevin dynamics makes the transition from work- 
ing correctly (at large /3) to failing (at small /3), a good test is to compare distributions of observ- 
ables computed at /i = using different initial conditions. In the absence of a chemical potential, 
if the imaginary parts of the complexified fields are initially zero, 1 = 0, the forces in the imagi- 
nary direction will be zero always, and therefore the configuration will remain constrained in the 
manner of real Langevin dynamics. This is called a cold start. Alternatively, the imaginary parts 
of the fields may be initialised randomly according to a Gaussian distribution, called a hot start. 

A good quantity for making such comparisons is the maximal force ^ max . When <j) 1 = 0, it 
will be constrained such that K™^ < 6/3. With complexified dynamics, there is no upper bound 
and the drift terms can fluctuate over several orders of magnitude JTT|] . In Fig. |3| distributions of 
the maximal force term for high and low /3 using hot and cold starts are plotted. 

In the absence of a chemical potential, the difference in behaviour between hot and cold starts 
is clear. At large /3 the complex Langevin (hot start) configurations are driven to the same equilib- 
rium distribution as the real Langevin (cold start) ones. The drift terms are bounded by 6/3 in both 
cases. However, at small /3, the drift terms from the hot start configurations fluctuate over several 
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Figure 3: Distributions of K m3X using hot and cold starts for = (top) and pi = 0.1 (bottom), lattice 
volume 8 3 , at two j3 values. Note the different vertical scales. 

orders of magnitude, whereas the drift terms from the cold start are bounded by 6/3. Since the cor- 
responding result for the action density is incorrect, this indicates that complex Langevin dynamics 
is not reaching a correct equilibrium distribution. At nonzero chemical potential, yi =0.1, the dis- 
tributions coincide because, even in the case of a cold start, the chemical potential forces the field 
into the complexified space. Consequently, this problem is not caused by the initial conditions, 
but by the dynamics itself At large /3, there are occasional occurrences of large drift terms 
(Fig. 3(d)| , inset), but the distribution is still dominated by the peak at 6/3. At small /3, this peak is 
absent and the maximal drift terms are substantially larger (Fig. |3(c)| ). 



7. Summary 

Complex Langevin dynamics offers the potential for a general solution to the sign problem. 
The XY model at finite chemical potential was studied using complex Langevin dynamics and 
compared to results obtained with the alternative world line formulation, which is free of the sign 
problem. It is necessary to integrate the stochastic equations using a dynamic stepsize, in order to 
eliminate the problem of runaway solutions that appear when using a fixed stepsize. The action 
density was computed in the region around \i 2 ~ 0. A comparison between results from 

using an imaginary chemical potential, yi 2 < 0, with those at \i 2 > shows that complex Langevin 
dynamics gives the correct result at large /3 but an incorrect result at small /3 . An analysis of the 
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relative difference between the action densities obtained using complex Langevin dynamics and 
the world line formalism over a wide region in the /3-jU plane shows that the failure of complex 
Langevin dynamics is strongly correlated with the phase of the theory. In the ordered phase at 
j6 > jS c complex Langevin dynamics works correctly, but it fails in the disordered phase j3 < p c . 

To investigate this difference further, we computed the distribution of the maximal force term 
K max at n = using two different initial conditions in which the imaginary parts of the fields 
were set randomly (hot) or to zero (cold). When complex Langevin dynamics yields the correct 
result (at large /3), the configuration from the hot start is driven to the same distribution as from 
the cold start. However, where it fails (at small /3), the distributions do not match. We conclude, 
therefore, that the failure of complex Langevin dynamics is not caused by the sign problem, but 
rather by an incorrect exploration of configuration space. This conclusion is supported by results 
that indicate that complex Langevin dynamics can evade the sign problem, see e.g. Refs. ^ 0, |L^ ]. 
An investigation into the incorrect exploration of complexified field space is currently in progress, 
building on the ideas put forward in Ref. JT3Q. 
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